Hard-thermal-loop resummed pressure of a degenerate 
quark-gluon plasma 



Rudolf Baier^ and Krzysztof Redlich^'^ 
^ Fakultdt fur Physik, Universitat Bielefeld, 
Postfach 10 01 31 , D-33501 Bielefeld, Germany 
'^Institute for Theoretical Physics, University of Wroclaw, PL-50204 Wroclaw, Poland 

(February 1, 2008) 

We compute the pressure of a finite density quark-gluon plasma at zero temperature to leading 
order in hard-thermal-loop perturbation theory, which includes the fermionic excitations and Landau 
damping. The result is compared with the weak-coupling expansion for finite positive chemical 
potential fj, through order and with a quasiparticle model with a mass depending on jj.. 
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In this letter we extend the recent work by Andersen, Braaten and Strickland by applying hard-thermal-loop 
O (HTL) perturbation theory at the leading order |^ to a relativistic quark-gluon plasma at finite and positive chemical 
potential and at zero temperature. These investigations are motivated by the observation of the poor convergence 

■ properties of the perturbative expansion of thermodynamic functions describing the quark-gluon plasma at high 
' temperature |^ and at high density 0,^. However, efficient resummation schemes have been recently developed in 

the case of scalar field theory [|| and of gauge theories like QCD which improve considerably the convergence 

■ of the free energy at high temperature. 

^ , In the calculation of the high temperature behavior of the pressure of a gluon gas ^ the inverse HTL-resummed 
gluon propagator is important. In the case of a degenerate quark-gluon gas at T = and non- vanishing chemical 
potential /i the improved quark propagator plays the central role. The ideal gas behaviour at T = and 7^ is 
determined by a quark loop 
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' where we state the result for a one-flavor (TV/ = 1) massless quark-gluon gas with TVc = 3 colors. The fermionic 
Matsubara frequencies are given by po = i{2n + l)TrT + fi. The T = limit is taken after the Matsubara sum is 
O ^' performed. 

D , The HTL-pressure Phtl receives contributions from three terms: Phtl — Pq + Pg + '^Pj that is from the HTL- 
■ resummed quark propagator, from the gluon propagator and from possible counterterms, which renormalize the 
pressure |jl| . The quark-loop contribution Pq is obtained by replacing the momentum in (|l]) by — Shtl , where 
. • Shtl is the quark-selfenergy in the HTL approximation . In the following we use the expressions and notations as 
[ given in i.e. 
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^ - Shtl = Am - Asjp with Aa=po- ^Oo(-) and ^ = po + ^(1 - -Qo(-)), (2) 

p p p p p 

where Qq denotes the Legendre function. The quark thermal mass at T = and /it > is given by 



M^ (3) 



with the MS renormalization scale jn. Therefore, for Nc = 3 we have to/ < ji. In terms of the functions Aq, As the 
pressure is 



Pq = 27V< 



^og{Al-Al), (4) 



where J]= ^X]n/^3 ''J (■iJ)'^ represents the sum-integral in d-dimcnsion, fi^ denotes the MS renormalization scale, 
necessary for the renormalization of the pressure using dimensional regularization for c? ^ 3 dimensions ^ . 

In order to evaluate (|4|), it is useful to use the spectral representation of the quark-propagator. Having in mind the 
linear dependence of Aq and Ag on to^ we first separate the ideal gas contribution by 
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For A±{po,p; m?) = (Aq T ^s) ^ the following dispersion relations hold 



A±{pq + iT],p;m ) = - — . (6) 



The sum rules for the HTL spectral function 



/ ^(^'o Tp)p±(po,p;m2) = 

turn out to be useful in our further discussion. Using (|^) and we find for the expression 
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Pq ^ Ptdeal + '2Nc I -z::2'l I iiPo ^ p) P+iPo , P] m^) + (po + p) P- {po, P] m^)} ■ (8) 



The derivation of is based on the observation that the expression in the curly brackets of (j^) is a meromorphic 
function in pq, which decreases faster than Pq^ for \pq\ — > oo, because of (|^). This allows to perform the Matsubara 
frequency sum and to take the T = limit in the standard way j9j. 

The representation covers the following contributions: (i) the fermion quasi-particle excitations for po > p, 
which are solutions of ReA^^(po = ^±,p) = 0. In (||) they contribute for uj±{p) > fi. (ii) the Landau damping 
contribution restricted hy fi < po < p. 

In order to proceed we state the explicit expressions for the spectral functions, 



p±ipo,P]m^) = 27r [Z±{p)dipo-io±ip)) + Z^{p)S{po+Lo^ip))] + -m\l T x)9{l - x'') 
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with X = PqIp^ L — log 1 and the residues at the quasi-particle poles Z±{p) = "^^^2 ^ ■ Next we consider the 
contributions (i) and (ii) separately and decompose 

Pq Pideal Pq, quasi ^" Pq.Ld ^" Pq.sub- (-^0) 

The quasi-particle contribution Pq^quasi to the pressure Pq is obtained as follows. Due to the (5-functions in the 
Po-integration in (||) is straightforward. We also note the identities dw± = (uj± Tp)Z±{p)^^, which allow to perform 
the TO^-integration. The quark mode contribution has an ultraviolet divergence for d=3 dimensions. Following 
the procedure of [0 this divergence is isolated by subtracting expressions from the integrands to obtain finite results 
for d — 3, and evaluating the subtracted integrals in c? = 3 — 2e, (e > 0) dimensions. The subtracted terms take into 

account that lj-i- p -\ — ^ — ^log^, for p » mj Q. The quark (— ) mode contribution is ultraviolet finite, 
therefore it does not require subtracted terms. The finite quasi-particle contribution for d = 3 reads 

f d^p 

Pq,quasi = 2A^c J J2nf^^^^^^^ ^ P')^{^+ ^ P')^{^^ - P) + (p) - p)0{uj- - fJ.)9{fl - p) 



+ i^-{p) -p)e{p-fi)- {u+{p) -p)e{fi~p) + 
2(p2 + 2m2)3/2 ^^^2 



+ + .,.2^ 42,3/2 a»g - m. (11) 



The Landau damping term Pq^Ld is treated in an analogous way. The dependence of p± on allows to perform the 
m^-integration in (||) explicitly. The ultraviolet divergence, which is also present in this contribution, is treated as 
described above. The finite term for d = 3 becomes 

Pq,Ld = -47Vc / ^(^^ - A^) ^ bo, p; m^) - $_ (po, p; m}) 

TTTJi'jPo 1 7rm J p + Po 

+ p3 (p2_p2 + 2m2) ^^^^ 
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where the angles ^± are given by 



7r(l =F x)!^^ 

$±(po,p;m.) = arctan ^ . (13) 

2p[p{lTx) + ^{2±ilTx)L)] 

The terms in ( |l^ ) which are subtracted in order to make it finite are constructed by observing that — — > 

4 

-^^[2a;/(l - + 2mJ/p2) + L] for p >> to/. 

Finally, we discuss the term Pq^sub, which contains the renormalized subtracted contributions evaluated in the limit 
e — > 0. We note that the quasi-particle and Landau damping terms in ^ have poles proportional to 1/e^ and to 1/e, 
which however cancel when adding these two contributions. Therefore no counter term is required: 



APg = 0. (14) 



The remaining finite term reads 



We observe from ( |15| ) that Pq includes terms of 0(my=a^). Consequently we have to consider contributions of 
0{mg=a^) arising from the HTL-resummed gluon loop. In leading order the result is already given in ||l[ when 
discussing the T ~ limit of the pressure, namely by 

Pq = ~{N^ - 1)-— TO4(log - 0.332837), (16) 

where the MS scale fis^ = 47re~'''^§ is introduced. However, here we have to replace the thermal gluon mass, which 
is TOg — ^i^t"^ in a hot gluon gas, by to^ — f^A^^ Q for the case of degenerate quark-gluon gas. We note that nig 
coincides with to/. In the following we fix the renormalization scale /I3 by /I3 = 0.717to/, such that Pg = 0, and only 
the quark-loop contributes to Phtl (besides the proper counter term for the gluon- loop). 

For later comparisons with the weak-coupling expansion we expand Pq in powers of ruf/fj, << 1, and we find 

P„ Tn?t frA n-o 

:^ 1 - 27V,^ + O(^) = 1 - 4^ + 0[al). (17) 

Pideal M M 

We note that the leading correction of 0(m^/^^) is exclusively due to the quark (+) mode contained in Pq quasi of 
®- 

The pressure of a degenerate quark-gluon plasma is calculated through order aj. in the weak-coupling expansion 
For one flavour {Nf = 1) it is given by 



°QCD — Pideal 



1 - 2 ' - (9.267 + log ( ^^^ ) log— )(^^)2 

TT TT 3 /Z4 TT 



(18) 



with £^ = ^(1 - |Mi2|i) for TV, = 3 and Z = log (mI/AH_). 

Although Pqcd in ((ll) depends on /i4 beginning at order a^, the dependence on turns out to be rather strong. 
For fi/Ajjg = 2, Pqcd / Pideal varies from 0.68 to 0.82 under the variation from ^4 — 2/i to 114 — 8/i. With increasing 
ji/ K-jjq the variation obviously decreases, and consequently the difference between the 0{as) and ©(a^) contributions 
in (llq) is diminished. This is shown in Fig.l for the choice 114 = 4/i. 

When comparing the 0(a<,) correction term of Pqcd in (|l8|) with the one of the HTL-resummed pressure Phtl 
at T = (^7|) we note a difference in magnitude by a factor 2. As in the case of the gluon gas at high T [Q HTL 
resummation overestimates the 0{as) term. In order to correct for this difference, we add a correction term to Phtl, 
namely Pcorr = 2^ Pideal- We define the pressure 

P = Phtl + Pcorr, (19) 

and attribute Pcorr to quasi-particlc interactions [Q. Presently, however, we are only able to approximate the corre- 
sponding series for Pcorr by its leading 0{as) contribution. 
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FIG. 1. The pressure for a degenerate quark-gluon plasma normalized to the ideal gas result as a function of yi/Kjjg. The 
Phtl + Pcorr pressure is shown for /X4 = 4/^ and as a band that corresponds to varying /i4 from /i4 = 2^ to ^4, = 8/^ (solid 
line). The weak-coupling expansions through order 0{as{li4,)) and 0(as(/i4)) are shown by long-dashed and dotted lines 
correspondingly; the HTL pressure Phtl is indicated by the dashed-dotted line, and the pressure Pmodei for non-interacting 
massive, m — m/(/i4), fermionic quasi-particles is shown as short-dashed line. All are calculated for ^4 = 4/i. 



In Fig.l we show Phtl as a function of /i/A-jg-j evaluated for ^4 = 4/i, indicating the overestimate of the 0{as) 
correction. Numerically the main contribution to Phtl is due to the quark mode to+ip) in ([Tl]), which is arising from 
the momentum interval fi > p > 0. The result for P is plotted in Fig.l for 114 = and as a band which corresponds to 
the variation of the scale by 2ii < ji^ < 8^. We observe that the O(a^) corrections of the weak-coupling expansion 
( |l8|) lead to a much bigger deviation from the ideal gas behaviour than the HTL-rcsummcd P, which follows more 
closely the 0{as) behaviour of ( p^ ) in the considered region of /j/Aj^g. 

We can not compare at present our result for the pressure at T = 0, /x > 0, with lattice results, as it could be 
done at /I = and T > 0. Instead we evaluate P for non-interacting massive, to ^ 0, fermionic quasi-particles 
{Nf = 1, Nc = 3) with finite chemical potential. The corresponding pressure reads ||l|] 
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(20) 



with /i > m. Choosing m = nif of (^ we realize that through 0{as) Pmodei coincides with Pqcd ( P^ ) as well as with 
P (19). This is clearly seen in Fig.l where Pmodei is plotted for /i4 = 4/x. Indeed the quasi-particle approximation 
^ turns out to be an excellent effective description of the degenerate quark-gluon plasma at T = 0, when the quasi- 
particle mass TO is identified with the quark thermal mass to/ of the HTL approach. This agreement is also shown 
in Fig. 2, where we summarize different approximations for the pressure discussed in this letter, which we normalize 
to the ideal gas for T = 0, to = and positive ji. For the choice /i4 = 4/i we plot in Fig. 2 the ratio P/Pideai as a 
function of the coupling as{n4) in the region < c^s < 0.4 for Phtl, P = Phtl + Pcorr, Pqcd through 0{as) and 
O(a^), and for Pmodei- 

Pisarski and Rischke in a recent paper JlSf conjecture that perturbation theory at T = 0, /i > may work better 
than for T > 0, i.e. they estimate its breakdown when the coupling becomes as large as 2.4. From Fig. 2 one 

may estimate a "critical" as ^ 0.4, i.e. where the term becomes as big as the as term of Pqcd- At P > a 
corresponding estimate gives as ~ 0.1 — 0.2. 

Only after the inclusion of, at least, the next-to-leading order correction in the HTL-perturbation scheme one 
could judge if the prediction for the pressure of a degenerate quark-gluon plasma becomes considerably more stable, 
especially when compared to the weak-coupling expansion (p^. Then it would be extremely interesting to compare 
with numerical lattice results for the T = finite density plasma which also, however, still have to be worked out. 
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FIG. 2. The pressure for a degenerate quark-gluon plasma normalized to the ideal gas result as a function of as{fJ,4). The 
contributions to the pressure correspond to those of Fig.l, consistently calculated with /14 — Afi. 

After having completed this letter we have found a relevant new paper by Andersen, Braaten, and Strickland [T^ ] 
where the quark contribution to the free energy of a hot (T ^ 0) but baryon free_(/i = 0) quark-gluon plasma is 
calculated to leading order in HTL-perturbation theory. In another recent paper pa the calculation of the baryon 
density in a self-consistent approximation in the sense of is presented as a function of fi/Aj^g. 
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